def aire_rectangles(f,a,b,n):
    longueur = (b - a) / n
    aire_inf = 0  
    aire_sup = 0 
    for k in range (n):
        aire_inf = aire_inf + longueur * f(a + k * longueur )
        aire_sup = aire_sup + longueur * f(a + (k+1) * longueur )
    return aire_inf , aire_sup
    
def f(x):
    return x**2


aire_inf , aire_sup = aire_rectangles( f, 0 , 1 , 1000 )
print("{} < Aire < {}".format(aire_inf,aire_sup))